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Abstract We extend the Abrams-Strogatz model for competition between 
two languages [Nature 424, 900 (2003)] to the case of n(> 2) competing states 

o ' 

c/3 ' (i.e., languages). Although the Abrams-Strogatz model for n = 2 can be in- 

terpreted as modeling either majority preference or minority aversion, the 
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two mechanisms are distinct when n > 3. We find that the condition for the 
coexistence of different states is independent of n under the pure majority- 
preference, whereas it depends on n under the pure minority aversion. We also 
show that the stable coexistence equilibrium and stable monopoly equilibria 
can be multistable under the minority aversion and not under the majority 
preference. Furthermore, we obtain the phase diagram of the model when the 
effects of the majority preference and minority aversion are mixed, under the 
condition that different states have the same attractiveness. We show that the 
multistability is a generic property of the model facilitated by large n. 

Keywords Consensus • Majority rule • Population dynamics • Social 
dynamics 

PACS 89.65.-s • 87.23.Gc • 89.65.Ef 
1 Introduction 

The consensus problem, in which we ask whether the unanimity of one among 
different competing states (e.g., opinions) is reached, and its mechanisms are 
of interest in various disciplines including political science, sociology, and sta- 
tistical physics. In models of consensus formation, it is usually assumed that 
each individual possesses one of the different states that can flip over time. 
The flip rate depends on the environment such as the number of the individ- 
uals that adopt a different state. Statistical physicists have approached this 
problem by analyzing a variety of models including the voter model, majority 



A model of competition among more than two languages 



3 



rule models, the bounded confidence model, Axelrod's model, and the naming 
game (see [3] for a review). 

A major mechanism that would lead to consensus in a population is pref- 
erence for the majority. Collective opinion formation under various majority 
voting rules has been examined for mean-field populations [51 [TTl[T^l[T51ll5lll{Jl 
|2"6] and in different types of networks such as regular lattices [71151151 1^1151] . 
small- world networks |21j . heterogeneous networks |17j . and networks with 
community structure 1161118] . The majority preference may be identified with 
the aversion to the minority. When there are only two states, they are equiv- 
alent because one state is the majority if and only if the other state is the 
minority. However, the two principles may be distinct when more than two 
states are assumed . We are concerned with this case in the present study. 

Language competition is an example of consensus problems. The model 
proposed by Abrams and Strogatz (AS model) accounts for extinction of en- 
dangered languages PQ. The AS model implements competition between two 
languages for speakers in a population. The dynamics of the model is based 
on the majority preference, which is also regarded as the minority aversion 
because there are just two competing languages. 

Several authors found that different languages can stably coexist in variants 
of the AS model. Two languages can coexist by spatial segregation in a model 
in which competition dynamics and spatial diffusion are combined |281I29] . A 
Lotka-Volterra variant of the AS model also leads to coexistence [3D] . Intro- 
duction of bilingual individuals also enables coexistence [2"Sll25| . Castello and 
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colleagues investigated variants of the AS model with bilingualism on various 
networks [4 , 5 , 32 , 33 . Coexistence also occurs when the attractiveness of lan- 
guages is dynamically manipulated [B] or when bilingualism, horizontal (i.e., 
from adults to adults) and vertical (i.e., from adults to children) transmission 
of languages, and dynamics of the languages' attractiveness are combined [2"3l 



In the present work, we extend the AS model to the case of competition 
among a general number of languages, denoted by n. Our model is a mean- 
field model (i.e., without spatial or network structure), as is the original AS 
model. Because the AS model has been used for modeling competition of other 
cultural or social traits such as religion [2], opinion, and service sectors [20] . 
we use the term "state" instead of "language" in the following. We show that 
the behavior of the model is essentially different between n = 2 and n > 3. 
In particular, the coexistence of different states and the consensus can be 
multistable, i.e., the coexistence and consensus equilibria are both stable, only 
when n > 3. 



We extend the AS model to the case of competition among n(> 2) states. The 
dynamics of the fraction of state i (1 < i < n) is given by 



[35] 



2 Model 




(1) 
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where Xi is the fraction of state i in the population, and Pji represents the 
transition rate from state j to state i. Equation (JTJ) respects the conservation 
law Xi — 1. The transition rates of the original AS model (i.e., n = 2) 

are given by 



where a > controls the strength of the frequency-dependent state transition, 
Si > is the attractiveness of state i, and s i — 1 [I]- Because c simply 

specifies the time scale of the dynamics, we set c = 1. 

Equation ^fy allows two interpretations: majority preference because Pji — 
SiXl and minority aversion because Pji = Si(l — Xj) a . The two principles lead 
to the same model when n — 2. However, the two principles are distinct when 
n > 3. Therefore, we redefine Pji to allow for independent manipulation of the 
two factors. The transition rates of the extended model are defined by 



where /3(> 0) and a — /?(> 0) represent the strength of the majority preference 
and the minority aversion, respectively. When n — 2, the dynamics given by 
the substitution of Eq. ((3]) in Eq. JT]) becomes independent of the f3 value. 

3 Analysis 

3.1 Case of majority preference (i.e., (3 — a) 

In this section, we set j3 = a to analyze the case in which the majority prefer- 
ence is present and the minority aversion is absent. Substitution of Eq. in 



Pji - csixf = c Si (l - Xj) a - (1, 2), (2, 1)) 



(2) 




(3) 
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Eq. (J) yields 

-j 71 n 

&Xi ^ > a ^ > _ 

~ x j s i x i ~ x i 2-^ S 3 X 3 

= {s^- 1 - (sx*- 1 )) x u (4) 

where (■) represents the average over the population, i.e., the average of a 
state-dependent variable with weight .r,- (1 < i < n). Equation (|3J is a repli- 
cator equation [TJ] in which SiX^ 1 and (sx a ^ 1 ) = Y^i=i s i x e P^ a y the ro ^ e °f 
the fitness for state i and the average fitness in the population, respectively. 
The dynamics given by Eq. Q has n trivial equilibria corresponding to the 
consensus, i.e., the monopoly of a single state, and an interior equilibrium 
given by 

i 

<= ^ _l_ (l<i<«). (5) 

V(x) = — (sx a_1 ), where x = (x±, . . . , x n ), is a Lyapunov function of the 
dynamics given by Eq. because 

dV(x) ^ a -ldXi 

i=l 
n 

i=\ 

= -a (((sx - 1 ) 2 ) - (sa; - 1 ) 2 ) < 0. (6) 

V has a unique global extremum at x* , which is minimum for a < 1 and 
maximum for a > 1 ( Appendix [A"| . Therefore, the coexistence equilibrium 
given by Eq. ((SJ is globally stable for a < 1 and unstable for a > 1. 

Equation ^ also admits a unique equilibrium for each subset of the n 
states. When n = 3, for example, the equilibrium in which states 1 and 2, 
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but not 3, coexist is given by Eq. ([5]) for i = 1 and 2, with the denominator 

i i 

replaced by sl~ a + S2~ a , and 2:3 = 0. In general, there are (", J equilibria 
containing n' states. If 2 < n' < n— 1, these equilibria are unstable. For a < 1, 
the instability immediately follows from the fact that x* is the unique global 
minimum of V{x). For a > 1, any equilibrium containing n' (2 < n' < n — 1) 
states is unstable because it realizes the global maximum of the same Lyapunov 
function restricted to the simplex spanned by the n' states. 

When a = 1, we obtain V(x) = — (s) = 5Z i=1 SiXi. Therefore, if Si > Sj 
U *)i the consensus of state i is eventually reached. If Si — > Sj (j ^ i, i'), 
for example, the n — 2 states corresponding to Sj (j 7^ are eventually 
eliminated. The dynamics then stops such that states i and i' coexist. If all 
the three Si values are equal, any population is neutrally stable. 

Figure Q] represents the dynamics in the two regimes with n = 3, which we 
obtained by numerically integrating Eq. For a < 1, a trajectory starting 
from anywhere in the interior of the phase space, i.e., x that satisfies x\ +X2 + 
X3 = 1, x±, X2, X3 > 0, asymptotically approaches the coexistence equilibrium 
(Fig.QJi). It should be noted that a point in the triangle in Fig.QJi corresponds 
to a configuration of the population, i.e., x. For example, corner (i — 1, 2, 
or 3) represents the consensus (i.e., Xi — 1 and Xj = (j =/= i)), and the 
normalized Euclidean distance from the point to the edge e2~e3 of the triangle 
is equal to the x\ value. For a > 1, a trajectory starting from the interior of 
the triangle converges to one of the n consensus equilibria, depending on the 
initial condition (Fig. [T}d). 
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In Fig. [21 a bifurcation diagram in which we plot x\ against a is shown for 
S\ = 0.40, S2 = 0.35, and S3 = 0.25. As a approaches unity from below, the 
stable coexistence equilibrium approaches the unstable consensus equilibrium 
corresponding to the largest Sj value (e± in Fig.[T]). At a — 1, the two equilibria 
collide, and an unstable coexistence equilibrium simultaneously bifurcates from 
the consensus equilibrium corresponding to the smallest Si value (e3 in Fig. [1} . 



Fig. 1 Dynamics of the extended AS model when the majority preference is present and 
the minority aversion is absent. We set n = 3, /3 = a, si = 0.40, S2 = 0.35, and S3 = 0.25. 
(a) a = 0.5 and (b) a = 1.4. Solid and open circles represent stable and unstable equilibria, 
respectively. 



3.2 Case of minority aversion (i.e., (3 = 0) 

In this section, we set /3 = to analyze the case in which the majority prefer- 
ence is absent and the minority aversion is present. Substitution of Eq. in 
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Fig. 2 Bifurcation diagram for the extended AS model when the majority preference is 
present and the minority aversion is absent. The parameter values are equal to those used 
in Fig. [l] except that we vary a. 

Eq. ((U yields 

j, — Si ^ ^ Xj ( 1 Xj ) ( 1 Xi ) 

^((l-xn-a-^r^. (7) 

In contrast to the case of the majority preference (Sect. 13. ip . the simplex 
spanned by n' (2 < n' < n — 1) states is not invariant under the dynamics 
given by Eq. (0. Therefore, a state that once gets extinct may reappear. In 
this section, we numerically analyze Eq. ([7]) for n — 3 and n = 4. For general 
n, we analytically examine the special case in which Si is independent of i in 
Sect. [331 

For n = 3, the dynamics for various values of a is shown in Fig. [3l We 
set si = 0.36, S2 — 0.33, and S3 = 0.31. When a is small (a < 1), there is 
a unique globally stable coexistence equilibrium in the interior (Fig. [3^,). The 
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three consensus equilibria e±, e%, and are unstable. At a = 1, e±, e2, and e% 
change the stability such that they are stable beyond a = 1 (Appendix [Bj . Si- 
multaneously, a saddle point bifurcates from each consensus equilibrium. The 
bifurcation occurs simultaneously for the three equilibria at a = 1 irrespective 
of the values of si, S2, and S3. Slightly beyond a = 1, the three consensus equi- 
libria and the interior coexistence equilibrium are multistable (see Fig. [3Jd for 
the results at a = 1.3). As a increases, the attractive basin of the coexistence 
equilibrium becomes small, and that of each consensus equilibrium becomes 
large (see Fig.[5J; for the results at a = 1.4). At a — a c \ « 1.43, the coexistence 
equilibrium that is stable for a < a c \ and the unstable interior equilibrium that 
bifurcates from at a = 1, where i corresponds to the largest Si value (i = 1 
in the present example), collide. This is a saddle- node bifurcation. 

Numerically obtained a c \ values are shown in Fig. |4^i for different values 
of Si, S2, and S3. A point in the triangle in the figure specifies the values of 
si, S2, and S3 under the constraint s\ + S2 + S3 = 1, s,- > (1 < i < 3). It 
seems that a c i is the largest when Si = 1/3 (1 < i < 3). Figure S] suggests that 
heterogeneity in Si makes a c i smaller and hence makes the stable coexistence 
of the three states difficult. When Si as 1 and Sj w (j ^ i), we obtain a c i a 1. 

When a is slightly larger than et c i, there are two saddle points in the 
interior. In this situation, one of the three consensus equilibria, which depends 
on the initial condition, is eventually reached (Fig. [3ji). However, the manner 
with which the triangular phase space is divided into the three attractive 
basins is qualitatively different from that in the case of the majority preference 
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(Fig. [TJd). In particular, in the present case of the minority aversion, even if X\ 
is initially equal to 0, the consensus of state 1 (i.e., ei) can be reached. This 
behavior never occurs in the case of the majority preference and less likely for 
a larger a value in the case of the minority aversion (Fig. [5^). 

The sizes of the attractive basins of the different equilibria are plotted 
against a in Fig.[SJi. Up to our numerical efforts with various initial conditions, 
we did not find limit cycles. A discrete jump in the basin size of the coexis- 
tence equilibrium is observed at a c \ ks 1.43, reminiscent of the saddle- node 
bifurcation. Interestingly, the attractive basin of the consensus equilibrium e± 
is the largest just beyond a c \. 

As a increases further, the second saddle-node bifurcation occurs at a = 
a C 2 ks 2.81, where an unstable node and a saddle point coappear (Fig. [3f). 
Logically, the sizes of the attractive basins could be discontinuous at a — a C 2 
because some initial conditions with small x\ might be attracted to e\ when a 
is slightly smaller than a C 2 and to or when a is slightly larger than a C 2- 
However, up to our numerical efforts, we did not observe the discontinuity, as 
implied by Fig. [5^i. 

Numerically obtained a C 2 values are shown in Fig.|4]3 for different values of 
si, S2, and S3. Heterogeneity in Sj makes a C 2 larger. In addition, a C 2 is equal to 
<z c i when si = S2 = S3 = 1/3. In this symmetric case, the three saddle points 
simultaneously collide with the stable star node at a = 1. Beyond a = 1, 
the equilibrium that is the stable star node when a < 1 loses its stability to 
become an unstable star node. The three saddle points move away from the 




Fig. 3 Dynamics of the extended AS model when the majority preference is absent and the 
minority aversion is present. We set n = 3, j3 = 0, si = 0.36, S2 = 0.33, and S3 = 0.31. (a) 
a = 0.9, (b) a = 1.3, (c) a = 1.4, (d) a = 1.5, (e) a = 2.6, (f) a = 2.9, and (g) a = 10.0. 
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Fig. 4 Dependence of (a) a c i and (b) a c 2 on si, S2, and S3 when the majority preference is 
absent and the minority aversion is present. A point in the triangle corresponds to a triplet 
(si, S2, S3), where si + S2 + S3 = 1 and Sj > (1 < i < 3). 

unstable star node as a increases. This transition can be interpreted as three 
simultaneously occurring transcritical bifurcations. 

The unstable node that emerges at a — a C 2 approaches x* — 1/3 (1 < 
i < 3) in the limit a —¥ 00, as shown in Appendix [Cj The three saddle points 
approach (xi,x 2 ,x 3 ) = (1/2, 1/2, 0), (1/2, 0, 1/2), and (0,1/2,1/2), as shown 
in Fig. [3g. This is a trivial consequence of the proof given in Appendix [C] 
Therefore, the heterogeneity in Si does not play the role in the limit a — > oo 
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such that the phase space is symmetrically divided into the three attractive 
basins corresponding to e\, e%, and e%. 

For n = 4, the relationship between a and the sizes of the attractive basins 
of the different equilibria is shown in Fig. [SJd. The results are qualitatively the 
same as those for n = 3 (Fig. [S^,). 

3.3 Symmetric case 

In the previous sections, we separately considered the effect of the majority 
preference (Sect. 13. ip and the minority aversion (Sect. I3.2p . In this section, 
we examine the extended AS model when both effects can be combined. To 
gain analytical insight into the model, we focus on the symmetric case Sj = s 
(I < i < n). Although normalization Sj = I leads to s< = l/n, we set 

s = 1 in this section to simplify the notation; s just controls the time scale of 
the dynamics. Then, Eqs. (fTJ) and (jSJ) are reduced to 

, n n 

Equation ([5]) implies that, regardless of the parameter values, there exist 
n trivial consensus equilibria and symmetric coexistence equilibria of n' (2 < 
n' < n) states given by x* — l/n', where i varies over the n' surviving states 
arbitrarily selected from the n states. 

Owing to the conservation law Y17=i x i ~ •"■> ^ ne dynamics are (n — 1)- 
dimensional. The eigenvalues of the Jacobian matrix of the dynamics at the 
coexistence equilibrium containing the n states are [n — l)-fold and given by 
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Fig. 5 Sizes of the attractive basins for different equilibria when the majority preference 
is absent and the minority aversion is present. The lines with legend "State i" represent 
the basin size for the consensus equilibrium of state i. The lines with legend "Coexistence" 
represent the basin size for the coexistence equilibrium, (a) n = 3, f} = 0, si = 0.36, 
s 2 = 0.33, and s 3 = 0.31. (b) n = 4, = 0, si = 0.28, s 2 = 0.26, s 3 = 0.24, and s 4 = 0.22. 
We obtain a c i ~ 1.43 and a c 2 ~ 2.81 in (a) and a c i ss 1.91 and a c 2 ~ 3.29 in (b). We 
calculate the sizes of the attractive basins as follows. First, we take the initial condition (xi, 
x 2, = (O.Oli, O.Olj, O.Olfc), where k > 1, i + j + k = 100, for n = 3 and (xi, ^3, 
x 4 ) = (0.05i, 0.05jr', 0.05fc, 0.05f), where i,j,k,l > 1, i+j + k+£ = 20, for n = 4. Second, we 
run the dynamics starting from each initial condition until the trajectory converges. Third, 
we count fraction of the initial conditions that converge to each stable equilibrium. 
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{nY ( l ~ nT 1 [(n - 2)/3 + a - n + 1], as shown in Appendix [Ql There- 
fore, the coexistence equilibrium is stable if and only if 



Similarly, we show in Appendix [B] that the consensus equilibria are stable if 
and only if 



Coexistence equilibria of n' (2 < n! < n — 1) states are always unstable 
(Appendix [Dj. 

Figure [5] is the phase diagram of the model in which the stable equilibria 
for given parameter values are indicated. The thin solid and dashed lines sep- 
arating two phases are given by Eqs. © and (ITU)) , respectively. A multistable 
parameter region exists when n > 3; Eq. (0) is reduced to a < 1 when n = 2. 
When n > 3, the multistablity occurs except in the case of the pure major- 
ity preference (i.e., (3 — a). The multistable parameter region enlarges as n 
increases. 

4 Discussion 

We analyzed an extended AS model with n states. We showed that the in- 
troduction of the minority aversion as compared to the majority preference 
changes the behavior of the model with n > 3 in two main aspects. First, differ- 
ent states stably coexist up to a larger a value with the minority aversion than 



(n - 2)f3 + a - n + 1 < 0. 



(9) 



a > 1. 



(10) 
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Fig. 6 Phase diagram of the extended AS model when all the Sj values are equal. The thin 
solid lines and the thin dashed line are given by Eqs. (0 and 1101 . respectively. 

with the majority preference. Nevertheless, it should be noted that a is the ex- 
ponent associated with different quantities in the two cases (Eq. ©)• Second, 
the multistability of the consensus equilibria and the coexistence equilibrium 
is facilitated by the minority aversion and opposed by the majority preference. 
We verified that the main results also hold true in the case of more general 
transition rates than Eq. ([3]), expressed as Pji = sj(l — Sj) 1 ~ 1 x^ (1 — Xj) a ~P 
(Appendix [Ej . 

Volovik and colleagues examined mean-field dynamics of a three-state opin- 
ion formation model with minority aversion |34j . Coexistence of at least two 
states occurs in their model even if a random choice term, equivalent to diffu- 
sive coupling, which the authors assumed, is turned off. This is because only 
the most minor state decreases in the number of individuals and the other two 
major states are equally strong in their model. In our model, the most major 
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and second major states have different strengths in attracting individuals, and 
the coexistence equilibrium is stable only for small a. 

Fu and Wang considered the combined effects of the majority preference 
and minority avoidance on coevolution of opinions and network structure |10| . 
They assumed that the majority preference is used for collective opinion forma- 
tion and the minority avoidance guides network formation. They showed that 
segregated groups, each composed of individuals with the same opinion, evolve 
when the minority avoidance is dominantly used. The coexistence of multiple 
states owing to the segregation has also been shown in spatially extended AS 
models [28 , 29 . We showed that coexistence of different states is facilitated by 
the minority aversion, not by the majority preference, and that stable coexis- 
tence occurs without segregation, or other spatial or network mechanisms. 

Nowak and colleagues analyzed a replicator-mutator equation as a model of 
language evolution [27]. They showed that coexistence of different grammars 
(i.e., states in our terminology) and a consensus-like configuration in which 
one grammar dominates the others are multistable when a learning parameter 
takes an intermediate value. If the learning is accurate, the consensus-like 
configuration becomes monostable. Our model and theirs differ in at least two 
aspects. First, the control parameter in our model is a, that is, the strength of 
the sum of the majority preference and minority aversion. Second, the stable 
coexistence in our model requires some minority aversion, whereas their model 
takes into account the majority preference but not the minority aversion. 
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A Unimodality of the Lyapunov function 

We denote the Hessian of V(x) by H(V), which is an (n — 1) X (n — 1) matrix owing to the 
conservation law X)2=i x i = 1- We obtain 



Hij{V) 



d 2 V(x) 
dxidxj 

o(l - a) (s z x1- 2 + SnzJT 2 ) (1 < i = j < n - 1), 
a(l - a)s n Xn~ 2 (1 < i ^ j < n - 1), 



(11) 



where x n = 1 — Y27=i x i an <l s „ = 1 — Y^=i s i- F° r anv x > H(V) is positive (negative) 
definite for a < 1 (a > 1) because 



n— 1 n— In— 1 

y T H(V)y = "(1 - a) Sl x1- 2 yl + £ £ a(l - a)s„<- 2 KW 

i — 1 i — 1 j — 1 

-1 /n-1 \ 2 

a(l -a) + s„<- 2 £ Wl 



(a SI), (12) 



for any (n — l)-dimensional nonzero vector y. Therefore, V(x) is strictly convex (concave) for 
a < 1 (a > 1) and hence has a unique global extremum at a:*, which realizes the minimum 
(maximum) of V(x). 
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B Eigenvalues of the Jacobian matrix in the consensus equilibria 

On the basis of Eqs. ifjQl and (0, the (n — 1) X (n — 1) Jacobian matrix (Jij) = (dxi/dxj) 
(1 < i, j < n — 1) is represented by 



and 



where 



lim 

X£ — tx^ 



Jij — lim (/jj li n kij + kin) 7^ j)) 

— >a;* 
x k = x k (kytj) 



(13) 



(14) 



(15) 
(16) 



and 



(17) 



By exploiting symmetry, we examine the Jacobian matrix at the consensus equilibrium 
jiven by Xi = (1 < i < n — 1) and x n = 1. In this case, Eqs. (1131) and (114ft are reduced to 



Ju = asi lim x a - 

:c;->0 1 



(18) 



and 



(19) 



respectively. Therefore, irrespective of the value of /3, the eigenvalues of the Jacobian matrix, 
denoted by A; (1 < i < n — 1), are given by 



Aj — t/jj — < 



-s„ (a > 1), 
Si - s„ (a = 1), 
+oo (a < 1). 



(20) 
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C Coexistence equilibrium when a — >• oo 



For general n, Eq. (7) implies that 

(1 Xj) X^ 



<(l-x)») (l<i<n) (21) 



(i-xj)»x; _ (i 



is satisfied at the equilibrium. In particular, we obtain 

(i^j)- (22) 
Equation J22j leads to 

lim f log(1 _ 2?) + ^ _ i2£fA = lim f log{1 _ x|) + !2£fl _ . (23) 

Because the third term on each side of Eq. H23H is negligible, we obtain 

lim l—^l ( fl V = l. (24) 

The coexistence implies that x| and x* do not tend to as a — » oo. Therefore, Eq. II 241 1 
leads to x* = x* in the limit a — ¥ oo. Therefore, x* = 1/n (1 < i < n). 



D Eigenvalues of the Jacobian matrix in the coexistence equilibria 
in the symmetric case 

On the basis of Eq. (SJ), the (n— 1) X (n — 1) Jacobian matrix (Jij) = (dx\/dxj ) is represented 
by Eqs. (1131) — ( 11711 with Sj = 1 (1 < i < n). At the coexistence equilibrium given by x* = 1/n 
(1 < i < n), Eqs. (1 13H and (1 14ft are reduced to 

J ii = n (fc-Z)= f^Y ^ l{n-2)j3 + a-n + l] (25) 

and 

= 0, (26) 

respectively. Therefore, the eigenvalues of the Jacobian matrix are (n — l)-fold degenerate 
and given by 

\i = Ju=(-Y(l-~ S ) 13 {(n-2)P + a-n + l] (1 < i < n - 1). (27) 
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By exploiting symmetry, we consider the coexistence equilibrium of n' states given by 
x* = (1 < i < n — n') and x* = 1/n' (n — n' + 1 < i < n). In this case, Eqs. JT3t and (fl4l 
with Sj = 1 (1 < i < n) are reduced to 



(28) 



/3 (1 - £) a ^ lim Xi ^o zf- 1 - (^) /3_1 (1 < i < n - n'), 

i 1 ~ ^Y'^ 1 IK - 2)/9 + a - »' + 1] (»- n'+ 1 < i < n- 1), 

(1 - ^r^ 1 K»' - 2)/3 + a - n' + 1] (n - n' + 1 < i < n - 1 
+ (^) ,? -/ 3 ^( 1 -^) a " ,3lim ^-0^- 1 andl<i<n-n'), ( 29 ) 
(otherwise). 
The eigenvalues of the Jacobian matrix are the diagonal entries because the Jacobian matrix 
is a triangular matrix. Therefore, we obtain 

(^) /3_1 09 >1). 

a, -./„ - { (!_ ^y- 1 ~i 03 = i), ( 3 °) 

+oo (/3 < 1), 



for 1 < i < n — n and 



A A. I — I (1- — 



[(«' - 2)j8 + o - n' + 1] 



(31) 



for n — n' + I < i < n — I. Equations (130 1) and (1 3 1 I t imply that at least one eigenvalue is 
positive unless a = /3 = 1. When a = f) = 1, all the eigenvalues are equal to zero. 



E Analysis of a generalized model 

Consider a variant of the extended AS model in which the transition rate given by Eq. J3]l 
is replaced by 

Pji = -7(1 - s j ) 1 -'xf (1 - xtf-P, (32) 
where < 7 < 1. Equation J3j is reproduced with 7 = 1. Here, the exponents 7 and 1 — 7 
represent the strength of the preference for a strong state and the aversion to a weak state, 
respectively. 
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In the case of the majority preference (i.e., /3 = a), Eq. Q with the transition rates 
given by Eq. (1321) has n trivial equilibria corresponding to the consensus and an interior 
equilibrium given by 



37. (33) 

e?=i (»r(i-»<) i - ,y ) , ~ 1 

There also exists a unique equilibrium composed of arbitrarily chosen n' states (2 < n' < 
n — 1). The Lyapunov function is given by 



((l-s) 1 ^) 



y(») = - 77T . U -yV < 34 ) 



y(a;) is a Lyapunov function because 
^ 1 1 dt 

1 = 1 

= a ({(s^/- 1 ) 2 ) ((1 - s) 1-7 ) - (s''x a ' 1 (l - s) 1 " 7 ) (s 7 ^- 1 )) (35) 



and 



d((l-8)i->) _ " o^d* 

1—1 

n 

= £(1 - s,) 1 " 7 [.7x? ((1 - s) 1 " 7 ) - (1 - s,) 1 " 7 ^ (sV 1 )] 
i=i 

= (sV^l - s) 1 " 7 ) ((1 - s) 1 - 7 ) - ((I - s) 2 ' 1 - 7 ^ (sV 1 ) (36) 



lead to 

dt {(1-S) 1 -T> tt ((l_ s )l-7) a + 1 

a ((A 2 ) (B) 2 - 2 (AB) {A) (B) + (A) 2 (B 2 )) 
~ (B) a+1 

a((A(B)-(A) B) 2 ) 

= " XT < 0, (37) 

(B) a+1 ~ 

where A = s J x a ~ 1 and B = (1 — s) 1-7 . 

We have no proof of the unimodality of V(x). However, we numerically verified the 

unimodality for n = 3. These results are the same as those obtained in the main text for 

7 = 1 rsect. nri) . 
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For any /3, the Jacobian matrix (Jij) = (dx\/dxj) is given by Eqs. d 1 3 II and 1141 . where 
kji = I3s'((l - s.Y-'x^Xjil - Xj) a -^ (38) 



l 3l = «7(1 - s 3 f-^(\ - x,Y-P-\\ -(a -(3 + l) Xj ]xf. 



(39) 



and 



x n = l-J2 x i- (40) 

i=l 

In the consensus equilibrium given by Xi = (1 < i < n — 1) and x„ = 1, Eqs. (1 1 3 II and 
d 1 411 are reduced to 



J U = as](l - s n ) 1 -~< Iim xf- 1 - 4(1 - s^) 1 " 7 
— J-0 



and 



(41) 



(42) 



respectively. Therefore, the eigenvalues of the Jacobian matrix are given by 

4(i -s^-y (a>l), 
A ' \s]{l _ s „)i-T_ 4(1 - Si ) 1-7 = 1), (43) 

+oo (a < 1), 

and the consensus equilibrium is stable if a > 1. 

In the case of the minority aversion (i.e., j3 = 0), when n = 3, the sizes of the attractive 
basins for different equilibria for 7 = and 7 = 0.5 are shown in Figs.[7ji and[7p, respectively. 
The results are qualitatively the same as those for 7=1 (Fig.[5ji). 
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